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Abstract 



> 

It is shown how a integrable mechanical system provides all the localized static solu- 
, tions of a deformation of the linear 0(iV)-sigma model in two space-time dimensions. 

The proof is based on the Hamilton-Jacobi separability of the mechanical analogue 
system that follows when time-independent field configurations are being considered. 
In particular, we describe the properties of the different kinds of kinks in such a way 
that a hierarchical structure of solitary wave manifolds emerges for distinct N. 

*S ' 1 Introduction 

> : 

We divide this Introduction into three parts: A. A brief history and the "state of the art" . 
B. New developments and results to be presented in this work . C. Scenarios of possible 
physical applications. 

A. 

Kinks are solitary (non-dispersive) waves arising in several one- dimensional physical sys- 
tems. Here, we shall focus on the relativistic theory of A-interacting scalar fields built on 
a space-time that is the (l+l)-dimensional Minkowski space R 1,1 . In this context, kinks 
are finite energy solutions of the Euler-Lagrange equations, such that the time-dependence 
is dictated by the Lorentz invariance: <j>K(x,t) = qx (775=1) ■ Thus, the search for kinks 
leads to the solving of a system of A-coupled non-linear ordinary differential equations and 
therefore becomes a very interesting problem in Mathematical Physics. 

The study of topological defects began as an area of research in field theory by the mid- 
seventies; see [0 . It was immediately recognized that defects of the kink type are in one-to- 
one correspondence with the separatrix trajectories between the bounded and unbounded 
motion of a mechanical system, for which the motion equations precisely form the non-linear 
system of differential equations mentioned above. The equivalent mechanical system is also 
Lagrangian and thus automatically integrable if N = 1. For N > 2, complete integrability 
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is generically non-guaranteed and the equivalence to a mechanical system is not useful. 
This circumstance has been emphasized by Rajaraman; see |2| pp. 23-24, and partially 
circumvented by himself: the trial orbit method allows one to guesstimate particular types 
of kink trajectories. 

There are, nevertheless, theories with N = 2-coupled scalar fields such that the equivalent 
dynamical system is completely integrable. The prototype of this kind of system is the 
MSTB model: in || this was proposed in the context of the search for non-topological 
solitons with stability provided by a U(l) internal symmetry. In Reference M, the model 
was considered as a classical continuum approximation to a ID crystal with a two-component 
order parameter and it was shown that the search for kinks in this system requires that a 
completely integrable dynamical system be addressed. Ito, in a seminal paper ||, showed 
the Hamilton- Jacobi separability of the system of non-linear differential equations. He found 
all the kink trajectories and explained a very peculiar kink energy sum rule. 

Very rich manifolds of kinks were discovered in two N = 2 field theoretical models, close 
relatives to the MSTB system, in a recent research performed by the authors of the present 
work 0. The investigation of kink properties in these models requires the analysis of the 
separatrix trajectories in two related dynamical systems which are type I and III respectively 
in the classification of Liouville bidimensional (N = 2) completely integrable systems, see 
i- 

In fact, on choosing between the four types of Liouville dynamical systems those that 
meet appropriate critical point structure, one builds an enormous list of related N = 2 field 
theoretical models exhibiting manifolds of kinks of growing complexity; see [KJ. The role 
of these models can be understood by noticing that the MSTB system is a deformation 
of the 0(2)-linear sigma model. Instead of spontaneous symmetry breaking of 0(2) by a 
degenerated S 1 vacuum manifold, the 0(2) symmetry group is explicitly broken to Z 2 x Z 2 by 
a mass term; only invariance under <p a — * (— l) Sab (j) a , for b = 1,2, survives. From the point of 
view of quantum field theory, this deformation is very natural because in (l+l)-dimensions 
infrared divergences forbid the existence of Goldstone bosons, according to a theorem of 
Coleman 0. Even if it is absent in the classical action, a mass term will be generated by 
quantum corrections. 

We interpret this as follows: in the parameter space of the N = 2 relativistic scalar field 
theories invariant under the Z 2 x Z 2 with generators mentioned above, and potential energy 
of the form 

there are at least two distinguished points. There is a choice of coupling constants such that 
there is explicit 0(2) symmetry, which is spontaneously broken. This is the linear 0(2)- 
sigma model. The other interesting point is the MSTB model where the explicit Z 2 x Z 2 
symmetry generated by (p a — > (— l) Sab (fia, for b = 1,2, breaks spontaneously to the Z 2 sub- 
group generated by 2 — > — </> 2 . The key observation is that the renormalization group flow 
induced by quantum corrections in the parameter space avoids the 0(2)-sigma system and 
instead leads to the MSTB model, which also offers a variety of kinks. All the other field 
theoretical models exhibiting an abundant supply of kinks also correspond to deformations 
of 0(2)-symmetric systems with potential energies that depend on higher powers of <p\ and 

02,1,1- 
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There are strong analogies with the Zamolodchikov c-theorem, [IC]: deformations in the 
space of (l+l)-dimensional field theories leading from conformal to integrable systems are 
the most interesting ones. We meet an analogous finite dimensional situation: replace the 
(infinite dimensional) conformal group by the 0(2) group and integrability of one system 
with infinite degrees of freedom by integrability of a bidimensional mechanical system. 

B. 

This paper is devoted to investigating the kink solitary waves of the deformation of the 
linear 0(iV)-sigma model that generalize the MSTB system to the case of iV-inter acting 
scalar fields. Non-linear waves in relativistic field theories with N > 3 scalar fields were 
sketchily described for the first time in Reference ]lT[. In this work, we offer a detailed 
analysis of this issue. The following points merit emphasis: 

• a) The dynamical system that encodes the solitary waves of the model as separatrix tra- 
jectories has N first integrals in involution and hence is completely integrable. Passing 
from Cartesian to Jacobi elliptic coordinates in the "internal" space, R , the dynam- 
ical system becomes Hamilton- Jacobi separable. All the kink trajectories, and hence 
all the solitary waves, are then found by a special choice of the separation constants. 

• b) Deep insight into the structure of the kink manifold is gained by focusing on the 
N = 3 case. There are three kinds of kinks: 1. A two-parameter family of topological 
kinks with three non-null components that are "generic", i.e. they are not fixed under 
the action of the Z 2 x Z 2 x Z 2 group generated by (p a — > (—l) 5ab (fia, for b = 1,2,3. 
2. Four one-parameter families of "enveloping" non-topological kinks, also with three 
non-null components. The four families are related through the action of one Z 2 x Z 2 
sub-group and, together, form the envelop of the separatrix trajectories. 3. All the 
solitary waves of the N = 2 MSTB model appear "embedded" twice; once in each 
plane containing the two ground states. Different Z 2 sub-groups leave these embedded 
kinks invariant. 

• c) The structure of the kink manifold of the O(N) system with both explicit and spon- 
taneous symmetry breaking repeats the patterns shown in the N = 2 (MSTB) model 
and its generalization for N = 3. There are also generic, enveloping and embedded 
kinks, although when N increases the complexity of the kink manifold also increases. 
For instance, the N — 1 kink manifold is embedded N — 1 times in the manifold of 
kinks of the deformed linear 0(iV)-sigma model. 

• d) In a remarkable system obtained from the generalized MSTB model by also allowing 
asymmetries in the quartic terms of the potential, only the embedded and enveloping 
topological kinks living on singular edges survive as solitary wave solutions. In this 



system, proposed in Reference [|l!| for the N = 2 case, the energy of all the above 
topological kinks is exactly the same. Together with vacuum degeneration, there is 
therefore kink degeneration, a phenomenon that deserves further analysis. 

C. 

Solitary waves of the kind that we are to describe play an important role in condensed 
matter physics. Phase transitions characterized through order parameters of the vector type 
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are understood in terms of the linear (or non linear) 0(iV)-sigma model. The order parameter 
is organized in the fundamental representation of O(N) and the system becomes non-linear 
when this iV- vector is forced to take its values in the coset space M. = 0(N)/0(N — 1). 
In (1 + l)-dimensional space-time, kinks are accompanied by the fermion fractionization 



phenomenon fll3[ ; this describes the continuous approximation to the bizarre behaviour of 
certain one-dimensional polymers such as poly-acetilene. When the spatial dimension is 3, 
as in the real world, kinks become domain walls which are thus related to theories involving 
spontaneous breaking of discrete symmetries. This happens in the hot Big Bang cosmology, 
where domain wall topological defects can be formed in a phase transition occurring in 



the expansion of the very early Universe; see |L4]]. More recently, domain walls have been 
characterized as BPS states of SUSY gluodynamics and the Wess-Zumino model, [pT5| . In all 
these cases there are sets of scalar fields, as in our system, that presents a variety of domain 
walls with different characteristics when seen from a 3-dimensional perspective. 

In quantum field theory, the linear 0(iV)-sigma model describes systems with sponta- 
neous symmetry breakdown to an 0(N — 1) sub-group and N — 1 Goldstone bosons in the 
particle spectrum. At the beginning of the sixties Gell-Mann and Levy analyzed low energy 
hadronic phenomenology by introducing an effective Lagrangian field theory of this type 



16| . Besides becoming the central element of current algebra, linear sigma models also enter 
fundamental physics in the Higgs sector of gauge theories for elementary particle physics, 
see report |17j for a comprehensive review (of the perturbative sub-sector). For instance, the 
linear 0(4)-sigma model corresponds to the Higgs sector of the electro- weak theory, while 
the 0(24) x 0(5) case provides the bosonic sector of the SU(5) Grand Unified Theory. 

Either considered on their own or forming part of Gauge theories, there are reasons to 
discuss deformations of the linear sigma model. In the phenomenological approach, pions 
are identified with the Goldstone bosons of the model; a deformation is then necessary to 
convert these massless excitations in pseudo- Goldstone particles, accounting for the pions 
light mass. Gauge theories are today found in the low energy limit of (fundamental) string 
theory. Even though deformations in the bosonic sector of gauge theories produced by small 
mass terms spoil renormalizability, the low energy features remain (almost) untouched and 
it is (almost) legitimate to trust them. 

Here, we shall search for domain walls when these mild deformations are performed in the 
linear 0(iV)-sigma model. It is precisely in this kind of model where the cosmological problem 
of wall domination is avoided [0. Moreover, the system has a rich manifold of topological 
and non-topological solitons, allowing for topological defects with "internal" structure and 
leading to the existence of defects inside defects, a situation that generalizes a proposal of 
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The organization of the paper is as follows: In Section §2 we discuss the particle spec- 
trum of the deformed linear 0(iV)-sigma model as well as the manifold of the solitary wave 
solutions of the system. Section §3 is devoted to the N = 3 case, which is described in 
full detail. We describe the situation of the generalized MSTB model for any N in Section 
§4 and briefly discuss the phenomenon of kink degeneration in the Bazeia system. Finally, 
some conclusions are drawn and some new prospects opened in Section §5. An appendix on 
elliptic Jacobi coordinates is also offered. 
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2 Kinks in the deformed linear 0(iV)-sigma model 



In a generic sense we understand "kinks" as the solitary waves of a relativistic (1+1)- 
dimensional scalar field theory. We shall stick to the standard definition of solitary waves; 
see § and §: 

A solitary wave is a non-singular solution of the non-linear coupled field equations of 
finite energy such that their energy density has a space-time dependence of the form: 

e(x, t) = e(x — vt) 

where v is some velocity vector. 

Given one iV-component scalar field, which is a map from the R 1,1 Minkowski space-time 
to H N , x(x, t) = t), Xi{ x i t), ■ ■ ■ , Xn(%, t)), the dynamics of the system is governed by 

the action: 

S = J d 2 y^~d,x-d»x-V{x)} 
Here, [i = 0, 1 are indices in the space-time and we shall use a = 1, 2, . . . , N to label 

N 

components of the field in the "internal" H N space in such a way that x ' X = 51 XaXa- 

a=l 

In R 1 ' 1 we choose the metric as q = \ \ ^ | and the Einstein convention will be used 

v - 1 J 

throughout the paper only for the indices in R 1 ' 1 . The potential energy density is: 

V{xi,---,Xn) = ^(x-X~ ^ +Ej^ 

where A, m and (3 a are coupling constants of inverse length. The linear 0(iV)-sigma model 
corresponds to the case f3 a = 0, Va, which exhibits maximum O(N) symmetry. We shall 
focus on the deformation of this system, which is maximally non-isotropic in the harmonic 
terms, i.e. (3 a ^ /3&, Va 7^ 6. Somehow, the deformation is natural from a quantum field 
theoretical vantage point as we shall explain later and, moreover, we shall stick to the range 
ft 2 < m 2 , Va, in the parameter space because in this regime the structure of the kink manifold 
is richer. 

Introducing non-dimensional variables x ~ * Vti ~ > anc ^ ^ ~* we fi 11 ^ °ur 

expression for the action to be: 



S = ^Jd 2 x S^d^-dy- Vfa, ...An)} 



1 2^1 

V(0 1 ,...,0^) = -(0-0-l) +^-^ 

a=l 

2.1 Configuration space and particle spectrum 

The Cauchy problem for the field equations 



dV 

D( ^ = -^r> a = l,2,...,N (2) 
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is fixed by choosing a "point" <p(x, to) £ Maps(R, R N ) in the configuration space C , and 

its "tangent", (j)(x,t ) G T^Maps(R, H N ), as initial conditions to solve the system (|J) of 
non- linear PDE . 

The configuration space itself is isomorphic to the space of finite energy static configura- 
tions; if <p(x, t) = q(x), C is the set of continuous maps q : R — > R^ (q = (qi, . . . , qjy)) such 
that the (static) energy is finite: 

thus, C = {<f(x)/i£[g] < +oc}. q(x) G C only if <f satisfies the asymptotic conditions: 

lim — - = , lim q a {x) = v a , Wa = 1, . . . , N (4) 

x— >±oo flj; ai— >±oo 

where v = (ui, . . . , Ujv) is a constant vector that belongs to the set M. of vectors annihilating 
V. We assume, without loss of generality, the following ordering in the space of parameters: 
(7i = < cr 2 < . . . < (Tjv < 1- -M is thus formed by two vectors 

M = {v ± = (±l,0,---,0)} (5) 

which are the absolute minima of V. 

We refer to M. as the vacuum manifold because in the quantum version of the theory 
points in M. are the expectation values of the quantum field operators a at the ground states 
( "vacua" ) of the system. The vacuum degeneration - i.e. the existence of more than one vector 
in M. - is related to the breaking of symmetry. Besides two-dimensional Poincare invariance, 
there is a "internal" symmetry with respect to the discrete group G = Z 2 x xZ 2 = N 
generated by (p a — > ( — l) Sab 4> a , for b = 1, 2, . . . , N, Va = 1, . . . , N. The vacuum manifold is 
the orbit of one element by the group action 

M = G/Hz± = Z 2 , = v± 

Htf± = Z 2 x ^Yr. 1 xZ 2 is the little group of the vacuum tr 1 . The generators of H^± are 
the transformations <fi a — > (— l) Sab (p a , for b = 1, 2, . . . , N, and a — 2, 3, . . . , N, so that 
survives as a symmetry group when quantizing around v^ 1 . We can understand the internal 
parity group G as the discrete "gauge" symmetry: in (l+l)-dimensions no dynamical degrees 
of freedom related to gauge potentials appear. 



Vectors in A4 are critical points of V satisfying 



-fr- - = and therefore constant 
solutions of the field equations (0). The plane wave expansion around ± (x, t) = 

k 

is a solution of @ if the dispersion relation 

S ab u 2 = 6 ab k 2 + M 2 ab (v±) , M 2 ab = J^-(^) (6) 
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holds. 

In the quantum theory, these plane waves become the fundamental quanta with mass 
matrix M^ fe (-y ± ) and one reads the particle spectrum at a chosen critical point of V from 
(Q). Because tP" e Ai are minima of V there are no negative eigenvalues of Mf b ( , u ± ) and 
the dependence on time of the plane waves around v^ 1 is bounded: e luJt . The choice of 
as the starting point of the quantization procedure "spontaneously" breaks the symmetry 
G = Z^ N of the action to H$± = Z2 , which is the remaining one that survives in the 
particle spectrum. 

In our model, we read the particle spectrum from: 

/ 4 ... \ 

4 



M 2 (v±) = — 



a 2 ... 



(7) 



V ... a% J 

Considering this system as a physical description of the continuum approximation to 
a one- dimensional crystal with an iV-component order parameter, the particle spectrum 
describes a single phase with N phonon branches. We see explicitly how the symmetry 
group G = Z% N is "broken" by the choice of the vacuum to the = Z2 ^ sub- 
group: the N phonon branches have different masses or "energy gaps". From the point of 
view of particle physics we can say that there are no tachyons; only a pseudo-Goldstone 
particle becomes a Goldstone boson if the corresponding a a goes to zero. 

It is interesting to see the model as a member of the family characterized by the potential 
energy densities: 

1 ( N \ 2 N a b 

V = 9 \ ^ab4>a4>b - 7 2 + Y -TT^afa 

£ \a,b=l J a,b=l A 

where a a b, <J a b and 7 2 are "bare" non-dimensional parameters. Ultraviolet divergences are 
controlled by normal ordering in the quantum theory, but the need arises to introduce a 
renormalization 'point' /i 2 , and the dependence of the renormalized parameters on fi 2 is 
determined by the renormalization group equation. One special solution, a specific renor- 
malization group flow, might lead to the "point" : 

a£0" 2 ) = S ab ; ^Gu 2 ) = 0; 7V) = 1 

in the space of quantum field theory models in the family. This point is the linear O(N)- 
sigma model which has G = O(N) as the (continuous) symmetry group. The vacuum orbit 
is, however, M. = 0(N)/0(N — 1) = S 1 ^^ 1 , the (N — l)-dimensional sphere, and thus there 
is no unbroken symmetry left: there are N — 1 massless particles. If the only modification of 
the renormalized parameters is to allow for non-zero values of cr^,(/i 2 ), a = b = r + 1, . . . , N 
there are still r — 1 Goldstone bosons. 

Coleman |J established that in (1 + l)-dimensions the infrared asymptotics of the two- 
point Green functions of a quantum scalar field forbids poles at uj 2 = k 2 \ there are no 
Goldstone bosons in (1 + l)-dimensions. It is thus impossible to reach the 0(iV)-sigma 
model or its deformation with the 0(r) symmetry spontaneously broken to 0(r — 1) in the 
renormalization group flow. The closest admissible points are the models characterized by: 

Q^ 2 ) = Sab, 7"V) = 1, = 0,a^b 
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a«{pL 2 ) = < a* (fx 2 ) =a 2 2 <...< a NN {ia 2 ) = a% < 1 

In this paper we shall focus on the case of maximal explicit symmetry breaking; i.e. 
when strict inequalities in the parameter space occur. Nevertheless, we shall comment on 
the allowed situation characterized by 

°i = < a l = ■ ■ ■ = < Vn+i = ... = <%,<...< °r k +i = ■ ■ ■ = <>n < 1 

when there is degeneration in the spectrum but no Goldstone bosons. Note that the genera- 
tors of the 0(r\ — 1) x 0(r 2 — T\) x ... x 0(iV — r^) symmetry sub-group are in the little group 
of the vacuum. The symmetry group is G = Z 2 x 0(r\ — 1) x 0(r 2 — r±) X . . . x 0(N — r^), 
H$± = 0(ri — 1) x 0(r 2 — r%) x . . . x 0(N — r k ) and the vacuum orbit \s M = G/H^± = Z 2 . 

2.2 Configuration space topology: kinks and dynamical systems 

The configuration space of the model is the union of topologically disconnected sectors: 
2 

C — |_| C a/3 ; thus, 7r (C) = Z 2 x Z 2 and |vr (C)| = 4 are respectively the zeroth-order 

a,/3=l 

homotopy group of C and its order. This comes from the asymptotic conditions @ and the 
continuity of the time evolution . There are topological charges defined for each configuration 
in C as: d I 

Qa = ~ dx ~T ~ = o 0a(+OO,t) - (j) a (-00,t)) 

2 J-00 ax 2 

It should be noted that Q T a is independent of t, Va, and in our system equal to zero if 
a > 2. Therefore the four sectors C af3 are labelled by the values a, ft of the fields at infinity 
compatible with finite energy and determines the homotopy class in ttq(C) = Z 2 x Z 2 . 

The critical points of E are time-independent finite-energy solutions of the field equations. 
If they are not spatially homogeneous, the critical points correspond to solitary waves that 
are therefore related to the topological structure of C Q/3 . Besides complying with @, solitary 
waves satisfy the system of ordinary differential equations: 

d 2 q a = dV_ 
dx 2 dq a 

Recall that <j) a (x,t) = q a {x). Solving the system (§) is tantamount to finding the solutions 
of the Lagrangian dynamical system in which x = r plays the role of time, the "particle" 
position is determined by q a (j), and the potential energy of the particle is U(q) = —V(q). 
From this perspective the static field energy E is seen as the particle action: 

Trajectories that behaves asymptotically in the r-time as ruled by (|J) have a finite action, 
J, in the mechanical problem and are in one-to-one correspondence with solitary waves/kinks 
that have energy E = J in the field theoretical system. 

The mechanical analogy is very helpful when one is dealing with a real scalar field theory 
because, then, a first integral is all that we need to find all the solutions. Vector scalar 
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fields of iV components lead to N- dimensional dynamical systems which are seldom solvable. 
Magyari and Thomas 0] realized that the two-dimensional dynamical system arising in 
connection with the MSTB model is a completely integrable one in the Liouville sense; there 
are two first integrals in involution. Moreover, Ito [H has shown that the mechanical system 
is Hamilton- Jacobi separable, finding all the trajectories and hence all the kinks of the MSTB 
model. In a recent publication |5]], we have developed this procedure for two N = 2 models 
with interesting features: the first system is a deformation of the (l+l)-dimensional scalar 
field theory, where the potential energy density is the Chern-Simons-Higgs potential arising 
in self-dual planar gauge theories. The second one is a deformation of the linear 0(2)-sigma 
model, which is different from the MSTB model. 

To extend this method of finding kinks to the linear 0(N)-sigma model, N > 3, deformed 
in such a way that the O(N) symmetry is explicitly broken to G = Z£ N , we start from the 
"particle" action: 



1 dq dq 1 t _ _ 1%2 , 1 

a I 

The particle motion equations are: 
d 2 q a 



dr 2 



2q a (q-q-l) + a 2 a q a , Va = l,...,iV (10) 



which are mathematically identical to the field equations for static configurations. Finite 
action trajectories, kinks in the field theory, should also satisfy the asymptotic conditions: 

lim = 0, lim q a (r) = ±5 al (11) 

T— >±00 Qff T^iOO 

We shall use the Hamiltonian formalism to integrate the mechanical system. The canon- 
ical momenta p a { T ) — §f- — ~^( r )i together with the positions q a (r), form a system of local 

coordinates in phase space. We should bear in mind that p a (r) = -$rr when going back to 
the field theory. The mechanical Hamiltonian 

11 N 1 

h = zP-P~2 ? - X ) 2 " E 2^ (12) 

a=l 

leads to the system of canonical equations 

dq a ( T dp a , , 

-T = {h,q a }, -r = {h,Pa} 
dr dr 

equivalent to (PD|). Given any two functions F(q,p), G(q,p) in phase space, the Poisson 
bracket is defined in the usual way: 

dF dG dF dG' 



N I dF 



Op a Op a dp ai 



Obviously = 0, but our mechanical system is full of other invariants. In fact, as early 
as 1919 Gamier solved the motion equations and described periodic trajectories in terms 
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of Theta functions: the kink trajectories of finite "action" correspond to a limiting case and 
are the separatrices between the periodic trajectories and unbounded motion. More recently 



Grosse, and other authors 211 have shown that the functions: 



N -y N 

b=l,byta a b °a 6=1 

Lb = p a qb - Pbq a 

are first integrals in involution: 

{h,K a } = {K a ,K b }=0 

There is a set of N + 1 invariants in involution: I\, K\, K2, ■ ■ ■ , K^. The dynamical 
system is not superintegrable, however, because there are only iV-independent invariants: 
K\ + K 2 + . . . + Kn = 2Ii + 1. According to the Liouville theorem, the iV-dimensional 
mechanical system is completely integrable and all trajectories can be found, at least in 
principle. 

At this point we pause to explain the singular nature of the deformation of the linear 
0(iV)-sigma model chosen from among many possibilities. The ^f4> 2 a terms explicitly break 
the 0(A)-symmetry of the linear sigma model; the case a\ = 0, Va = 1,2, . . . , N. In 
the mechanical system the O(N) internal transformations become ordinary rotations. The 
angular momentum components, l a b, conserved in the limit a\ = 0, Va, are no longer 'time'- 
independent if <j\ 7^ 0. There are, however, N invariants K a , which in the limit a a = 0, Va, 
are given in terms of the 0(iV)-invariants: the r Casimir invariants and the r generators of 
the Cartan sub-algebra, where r = | or if iV-even or -odd is the rank of the group. 
A warning: in the N = odd case, the energy must be added to the other N — 1 invariants 
built from the Cartan sub-algebra and the Casimir invariants. For any N, the maximally 
asymmetric chosen deformation is special because it retains enough symmetry to solve the 
mechanical system. There is no Lie algebra associated with K a however; since the invariants 
are quadratic in q a , p a , the action of K a in the phase space, given by {K a , qb} and {K a ,pb}, 
is non-linear. 

In (l+l)-dimensional field theory, the energy-momentum tensor: 

dC 

is divergenceless due to invariance under space-time translations. P M = / dxT 0>1 are thus 
conserved quantities whatever the values of a\. The O(N) "isospin" currents however, 

N 

Ja = Yl C abc<Pbd^<p c 
6,c=l 

are only divergenceless if a a = 0, Va. The c a bc are the Lie O(N) structure constants and 
the charges Q a = J dxJ® are not conserved if there is no symmetry with respect to the 
transformation generated by them. For static configurations, we have 

dxT 00 = E = J, T 10 = T 01 = 0, T 11 = h 
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N 



Ja — 0) Jl — E C a,bdbc 



>,c=l 



In terms of the 'isospin'currents, the invariants K a can be written as: 

k«= E -^■(e^ 1 ) fE^) + (^) 2 + (2--M-^E' 

6=l,6^a a b ~~ a a \c=l / \d=l ) \ OX J 6=1 



We expect that the time-evolution occurs in such a way that there is some equation of 
non-linear character 

F (dLadKa\ =Q 

\ dt ' dx J 

between K a and 

La= E -^■(f:^) fEc^°) + f^) 2 + (2- ( x a 2 )^-^E0 fc 2 

6=l,6^a °a \c=l ! \d=\ ! \ 01 ) 6=1 

which reduces to ^f- = ^ when a a = 0, Va. The situation is analogous to that occurring 
between conformal field theories and models with infinite-dimensional algebraic symmetry 
as in (l+l)-dimensional Toda field theories and Toda affine models [2^|. There are two 
differences: (1) the conformal group is infinite dimensional in (l+l)-dimensions. We have 
only one finite-dimensional group O(N) and thus we can solve only the static limit of the 
field theory model. (2) Due to the non-linear character of the deformation of the O(N) Lie 
generators, we do not even have a finite-dimensional Lie algebra. 



2.3 The Hamilton-Jacobi equation and kink trajectories 

The K a invariants defined in (^) are quadratic in the momenta, but not orthogonal (they 
contain terms in p a pb,a ^ b). Therefore, the Stackel theorem can not be applied to assure 
Hamilton-Jacobi separability. This problem is surpassed in our dynamical system with the 
choice of some suitable system of coordinates. The appropriate system is provided by elliptic 
Jacobi coordinates, with a choice of separation constants determined by the deformation 
parameters giving mass to the Goldstone bosons; a 2 a = 1 — o\ , Va = 1, 2, . . . , N. Thus we 
define: 

N 



a 2 = — — = " v " a/ (14) 

y a at Ml-1\ \ > 



A(^) 



N 



n 



b=l,b^a 



ruling the change of coordinates from Cartesian, q = (g 1; . . . , q N ) to elliptic A = (Ai, . . . , A^v)- 
In the Appendix, it is explained how the elliptic variables are split: 

- oo < Ai < a 2 N < A 2 < < . . . < a\ < \ N < 1 (15) 



Notice that formula (fT4|) coincides with formula (^) in the Appendix if we change q a by 



qN-a+i and choose r N ^ a+1 = ct 
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Together with formula flHD, this splitting means that the change of coordinates produces 
a map from a sub-space of R^, characterized as the set of points which are not invariants 
under the Z% N group generated by q a — > (— l) Sab q a ,b = 1, •••,JV, to the interior of the 
infinite parallelepiped Pjy (oo) obtained by replacing the inequalities in ( |i5|) by equalities: 
— oo <Ai<ct^<...<ct|<Ajv<1- Notice that in this map 2 N regular points in H N go 
to a single point in the interior of Pn(oo); Singular points lie in the R m , m = 0, 1, . . . , iV— 1, 
sub-spaces that are invariant under the action of some non-trivial element of G = Z xN . 
These singular sub-spaces are mapped into the boundary of Pjv(oo). 

The standard length of an interval in Euclidean space is expressed in elliptic coordinates 
in the form 

N N 

ds 2 = dq a dq a = ^ g aa (\)d\ a d\ a 

a=l a=l 

because the metric g aa (X), as derived in the Appendix, is: 

9aa{X) = ~\j(KY 9ab{X) = °' Va ^ b 

N 

where A(X a ) = JJ(A a - and Q 

6=1 



N 

/a(A) = / a (A 1 ,A 2 ,--- ! A JV )= n(Aa-A.) 



6=1 



Therefore, the Lagrangian reads: 



L = -U(qu...,q N ) 



2 a t! V>, 

i N / j\ \2 

= ^E^M^r) - ^a(Ai, . . . , Aat) (16) 
where the potential in elliptic coordinates is: 

n(\ f l Ar-( a -l)A: + (l-a + /?)Ar 

^a(Ai, . . . , A N ) = - 2^ 7) T7T\ ^ > 

a=l Z Ja(Aj 

N N N 

a = S ^ ' ^ = 2 2 ^ 

a=l a=l,a<6 6=2 

The computation of £/a(A) is highly non-trivial and requires the use of formulas that 
follow the Jacobi Lemma, such as ([76]), ([77|), (|T5|), etc. 

AT 

1 The standard notation in the literature on elliptic Jacobi coordinates is A'(A a ) = ]^[(^a — ^b), see 

6=1 

Appendix. We shall use / a (A) in the main text instead of A'(A ), to stress the fact that this quantity 
depends on all the components of A 
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The canonical momenta associated to the A„ variables are: 



9ab[X) — = 0aa(A) — 



6=1 

and, through the standard Legendre transformation, we write the Hamiltonian: 

Z 0=1 /a (A J 

The key point is that H can be written in Stackel's form: 

N TT 

H = £_fl£- = (19) 

a=l fa(X) 

» [-2A(K)*l ~ \ (A^ +1 - (a - 1)A? + (1 - a + P)X^)' 

a=l fa (A) 

such that the Hamilton- Jacobi equation 

is completely separable. We now prove this last statement. 

Fixing H = Ii, the first integral of energy, and having in mind the expression fll9|) of H, 
we write the solution of (|20| ) as: 

AT 

5 = -/ 1 r + ^5 a (A a ). (21) 

0=1 

Therefore, (pO|) reduces to : 

The Hamilton- Jacobi PDE equation ( pOl) becomes equivalent to the system of non-coupled 
ordinary differential equations 

H a ( X a ) = VX~' + r, 2 \ N a ~ 2 + ■■■+ VN-lXa + VN (23) 



d\ a 



where 



H a (g, A.) = -2A(X a ) (g) - \ (A^ 1 - (a - 1)A? + (1 - a + ^A^ 



(24) 



due to the identity 



AT \JV-1 ^ \N-2 N i 

a = /i ^ + % x: + • • • + vn x: -==*- (25) 

a=l /a(A) a =l /a(A) a=l /a(A) 
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-observe that 771 = I\- which follows from the Jacobi Lemma and the subsequent relations 
(Appendix) 

N \N-1 N \N-i 

^^ = 1, £^r- = 0,Vi = 2,...,iV 

a=l fa (A) o=l fa (A) 

Alternatively, one could take a more direct approach to show formula fl23"l) . We start from 
fl22|), written explicitly as 



#i(Ai 



(Ai-A 2 )(A 1 -A 3 )---(A 1 -A JV ) 



+ 



+ 



H 2 (X 2 



+ 



(A2 — Ai)(A2 — A3) • ■ • (A2 — Aat) 

Hn(^n) 

(Xn — Ai)(Ajv — A2) • • • (Aat — Aat_i) 



h (26) 



Here, each H a (X a ) is of the form ( p4[ ) due to the ansatz (|2T|). Multiplying (|26| ) by (Ai — A2) 
and setting Ai = A2 one sees that -Hi (A) = H 2 (X), and hence by symmetry, all H a (X), 
a = 1, 2, - • • , N are identical. It suffices therefore to find i?i(A). Multiplying ( pB[) by /i(A) 
one obtains 

i?i(Ai) + P 2 (A)^2(A 2 ) + • • • + P N (X)H N (X N ) = I 1 (X 1 - A 2 )(A X - A 3 ) • • • (A x - A*), (27) 

where each P a (X), a = 2, 3, • • ■ , N, is a polynomial of degree AT — 1 in Ai. Differentiating 
(|27| ) A" times with respect to Ai yields 



d N H 1 (X 1 ) 



dX{ 



N 



mix 



(2? 



whence it follows that Hx(X) is a degree A" polynomial in A with leading coefficient Ji, in 
agreement with equation (p3|) . 

The set of separation constants 77^, z = 1, 2, . . . , N is another system of first integrals in 
involution. They can be expressed in the elliptic phase space T*P/v(oo) as functions of A a 
and 7r a = ^ by solving the linear system of equations (p3|) in the unknown 77^ which is a 
Vandermonde system. The N — 1 roots of the polynomial 

X^ 1 + fx^ 1 + ... + f = (Xa-F 2 )(X a -F 3 )...(X a -F N ) 

together with the energy F\ = Ii form another system of invariants in involution. Both 
systems are related through the identities 771 = I\ and: 

Va = (-l) a Ix E F h F h ■ ■ ■ K, ia = 2, 3, . . . , N 

il<i2<—<ia 

Therefore, all the separation constants r) a are proportional to the "particle" energy Ix- 
Defining the polynomial B(X a ) in the form: 



B(X a ) = Af +1 - (a - 1)X N + (1 - a + (3 + 2r ]l )X^ 1 + 2 V2 X^ 2 + ■ ■ ■ + ^ N 
the solution of the differential equation (p3|) is a quadrature: 

1 . ( dSa 



Sa(X a ) 



-sign 



dX a J J \ 



B(X a ) 
A(X a ) 



dX r , 



(29) 
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and the general solution of the Hamilton- Jacobi equation reads: 



s= -" iT+ s? ign O/\ 





B(X a ) 




A(X a ) 



dX r , 



(30) 



The explicit integration of the quadratures in (|29|) requires the theory of Theta functions 
of genus depending on N. The action of the associated trajectories is infinite because they 
are either periodic or unbounded. The asymptotic conditions (f|) that guarantee finite action 
to continuous trajectories satisfying them also require that the energy used by the particle 
in these trajectories should be zero. This is so because I\\ T=±00 = and, being an invariant 
of the evolution, Ii = 0, Vr. 

We recall that the trajectories of finite action in an evolution lasting an infinite time are 
the kinks of the parent field theory system: one just trades the finite action of the trajectory 
for finite energy of the non-linear wave. Therefore, the kinks are the trajectories obtained 
when all the separation constants in ( |23| ) are zero: r] a = 0,Va. These are the separatrices 
between bounded and unbounded motion and the integrals in (|29"D are easier to compute. 

The explicit trajectories are also provided by the Hamilton- Jacobi principle, through the 
set of equations: 

dS 

la = t— , a = 1,2,..., N 
or) a 

where the 7 a are integration constants. In the hypersurface of the phase space determined 
by 771 = • • • = TjN = 0, the first equation 



N 



71 = - r + 9 si g n K) 

a=l A 



X a dX a 



A{X a ) 



a 



_l) A AT + (l_ a + W -l 



(31) 



rules the time-dependence of the particle in its journey through the orbit. From the field 
theoretical point of view, it provides the kink form factor. The other N — 1 equations, 
z = 2,3,...,N, 



El ■ / n f l dX a 



A(X a ) 



X^-(a-l)X» + (l-a + (3)X»- 



(32) 



determine the orbit in P/v(oo), the intersection of iV — 1 hypersurfaces in the configuration 
space. Therefore, there is a N — 1-dimensional family of kinks parametrized by the finite 
values of 7,. 

Although ( |3lD and (^) identify all the separatrix trajectories of the mechanical system 
and henceforth all the kink solutions of the deformed linear 0(iV)-sigma model, an explicit 
description of such solitary waves is difficult for two reasons: (1). (|3lD and (]32|) form a 
system of transcendent equations of impossible analytical resolution. (2). Even if it were 
possible, expressing back the solution in Cartesian coordinates through (|72|) for iV > 3 is 
another impossible task by analytical means. 
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3 N=3 



To gain insight into the nature of the different kinks of the model, in this Section we shall 
address in full detail the N = 3 case. We shall deal with a (l+l)-dimensional field theory 
including three scalar fields which transform according to a vector representation of the 0(3) 
group. The structure of the solitary wave solutions of the N = 3 system is extremely rich 
from different points of view and shows the behavioural pattern of the general case with 
iV-component fields. 

3.1 The general solution of the Hamilton- Jacob i equation 

The Hamiltonian of the underlying dynamical system reads: 



H 



( 2 , 2 i 2 \ 1 / 2 i 2 i 2 i\ 2 _ 

2 gW + fc + 93-1) - TT ( /2 



2 „2 



a 



3 „2 



(33) 



2 1A 2 

in Cartesian coordinates. To write the Hamiltonian in elliptic coordinates, note that for 
iV = 3 we have: 

a = l + aj + al , (3 = ajal + o\ + ct 2 
A(A a ) = (A a -l)(A -<7 2 )(A a -<7 3 2 ) 
A'(Ai) = (Ax - A 2 )(Ai - A 3 ); A'(A 2 ) = (A 2 - A X )(A 2 - A3); A'(A 3 ) = (A 3 - AO(A 3 - A 2 ) 
3 1 

Hence, H = ^ ——H a , where 

0=1 A {*a) 



-2A(\ a ] 



7i„ 



^l(^a - 02) (A a - O- 



(34) 



The separatrix trajectories, those in one-to-one correspondence with solitary waves of 
kink type in the encompassing field theory, are fully determined by the equations (|32"D re- 
stricted to the N = 3 case: 



Co 



\/l - Ai + o 2 



\/l - Ai - CT 2 



Vl — A 2 + cr 2 



vr^ 



0-2 



yl — A3 + a 2 



\/l - A 3 - a 2 



o-3sign(7ri) 



o-3sign(7r 2 ) 



<T3sign(7r3) 



\/l - Ai - cr 3 



VI - Ai + cx 3 



\/l - A 2 - 0-3 



\/l — A 2 + cr 3 



\/l - A 3 - 0-3 



\/l — A 3 + 0-3 



<T2sign(7ri) 



<T2sign(7r 2 ) 



o-2sign(7T3) 



where C 2 = exp{27 2 cr 2 (X3(cr2 — °" 3 )} is constant, and: 



(35) 



v/T^AT- 1 



0-20-3 (cr| - cr| )sign(7Tl ) 



VT^ 

VI - Ai + 0-3 



1 

0-3 



VT Tr AT + 



<72 



VI — Ai — cr 2 



<T3cr|sign(7ri) 



CT20-^sign(7Tl) 



yr^A,- 1 
vr^A2~+ 1 



0-20-3 (oj - o| )sign(7r 2 ) 
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a/1 — A~2 + 0"2 





-A 2 - 




Vi 


-A 3 - 


1 


Vi 


- A 3 + 


1 


Vi 


-A 3 - 


0"3 



CT 3 o-|sign(7r 2 ) 



0"3 



y/1 - A 2 + 0"3 



o-2CT 3 (o-|-o-|)sign(7r 3 ) 



CT20-|sign(7r 3 ) 



\/l — A3 + cr 2 



<72 



(36) 



Vi - A 3 + 0-3 

with C3 = exp{27 3 a2a3<7 2 ! af (erf — 

Integration of (|3l|) in the iV = 3 case shows the time-table of the particle in each trajec- 
tory, or, the kink form factor: 



Ci(r) 



Vi 


-Xi 


- cr 2 


o- 3 CT|sign(7ri) 


Vi 


-A! 


+ ^3 


Vi 


-Ax 


+ cr 2 


<r 3 o^ sign (712) 


Vi 


-Ai 


- cr 3 


Vi 


-A 2 


+ cr 2 


Vi 


-A 2 


- 0"3 


Vi 


-A 2 


- cr 2 


cr 3 o-|sign(7r 3 ) 


Vi 


-A 2 




Vi 


-A3 


- cr 2 


Vi 


-A 3 


+ ^3 


Vi 


-A3 


+ cr 2 




Vi 


-A 3 


- 0-3 



o-2cr|sign(7r2) 



o-2cr|sign(7r 3 ) 



(37) 



if Ci(r) = exp{2(7! + t)(o| — o"|)cr 2 0- 3 }. Therefore, there is a family of kinks parametrized by 
the integration constants 72, 73: it corresponds to the family of curves in P 3 (oo) determined 
by the intersection of the surfaces defined by ( |3"oD and (PE|). The third constant, 71, fixes the 
center of the kink, the point where the energy density reaches its maximum value. 

Better intuition of the kink shapes requires an interpretation of the solutions described 
by equations ( p5|) and (|36"D in Cartesian coordinates. We shall describe how is this achieved 
in the next sub-sections, but before this it is convenient to note some details of the change 
of coordinates from Cartesian to elliptic in R 3 : 



4 



4 



J_ ( l_ Al) (l_ A2) (l_ A3 ) 
02°3 



o\[al-al) 
-1 



^-A 1 )(^-A 2 )(^-A 3 ) 



(38) 



a 3 2 -A 1 )(^-A 2 )(^-A 3 ) 



'-2 



-2 



(To 



The change of coordinates is singular at the three R 2 coordinate planes; q\ = 0, g 2 = and 
g 3 = 0. The image of the q\ = plane is a unique face, A3 = 1, of the P 3 (oo) parallelepiped: 



- 00 < Ai < cr 2 < A 2 < a\ < A 3 < 1 



^2 

2) 



The g 2 = plane, however, is mapped into faces A 2 = and A3 = cr 
g 3 = plane goes to faces A 2 = erf and Ai = o\ of ^3(00 



(39) 
while the 



022(Al,Cr 3 ,A 3 ) = 022(Al,<72,A 3 ) = ^33 (Ai, A 2 ,cr^ 



Observe that <?ii(<t 3 , A 2 , A3) = 
#33^1, A 2 , 1) = 00. The whole R 3 space 
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is mapped in P 3 (oo). Due to the symmetry under the group G = Z 2 x Z 2 x Z 2 generated 
by la — » —Q a , the mapping (|38"D is eight to one in regular points of R 3 : to any point in the 
interior of P3(oo) correspond eight points in R 3 away from the coordinate planes. These 
planes are fixed loci of some subgroup of G. 

The asymptotic conditions @ in q a restrict the motion to the compact sub-space D 3 of 
R 3 bounded by the tri-axial ellipsoid: 

?? + ff + 4=l (40) 

and are satisfied by finite action and zero energy trajectories. Elliptic coordinates are best 
suited for demonstrating such a restriction. In this coordinate system D 3 is mapped to the 
finite parallelepiped Ps(0): 

< Ai < a\ < A 2 < o\ < A 3 < 1 (41) 

The unique non-singular face of ^(0) with respect to the change of coordinates is Ai = 
and the inverse image of this face is the ellipsoid (|40|) . The asymptotic conditions ([|) force 
r/b = 0,Vfe, and thus, by Q2lf), H a = 0,Va, for the finite action solutions. A 2 and A3 are 
bounded, see (|4~T|). Thus, we focus on, 

Equation (|42"D describes the motion of a particle with zero energy moving under the influence 
of a potential 



2 



Ai 



-2 



V(Ai) = 4^3T> -00 < A x < a 3 

= 00, a\ < Ai < 00 

V(Ai) has a maximum at Ai = and goes to —00 when Ai tends to —00; therefore, bounded 
motion occurs only in the Ai G [0, a|] interval and the trajectories giving rise to kinks lie in 
P3(0), seen in elliptic coordinates, or D 3 in Cartesian space. 

In Figure 1 the whole picture is depicted and we notice the following important elements 
of the dynamics: 

- Points: (1) the origin. This is a fixed point of G = Z x3 and thus only one point O in D 3 is 
mapped to the vertex O in P3(oo). (2) Points B, C, D: these are the intersection points of 
the three distinguished ellipses , g? + % = 1, g? + % = 1 and % + % — 1, in the ellipsoid 

(|40|). They are fixed points under the action of a sub-group Z 2 2 of G and thus, two points 
in the boundary of D 3 are mapped to a single point in the boundary of Ps(0). D is the point 
where the two vacuum points tP are mapped and hence it is a very important point of the 
dynamics: every finite action trajectory starts and ends at D. (3) Points Fx, F 2 , F 3 , the foci 
of the above ellipses. Again two points in D 3 are mapped in a unique point in Ps(0). (4) The 
umbilicus A of the ellipsoid Ai = is another characteristic point; in Cartesian coordinates 
A corresponds to four points in the boundary of D 3 because they are invariant only under a 
Z 2 sub-group of G. 

- Curves: 
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Figure 1: (a) The domain D 3 in R 3 : Cartesian coordinates (b) The domain Ps(0) in R 3 : Jacobi 
elliptic coordinates 



The ellipse with foci F 2 



4 + ^S = l (43) 

in the qs = plane passing through F 3 and Fj. This is the edge A 2 = cxf = Ai in -Ps(O). 
Observe that four points on the ellipse (|43|) are mapped to one point in the edge of 
Ps(0), because it is invariant under a Z 2 sub-group of G. The map leading to F 3 and 
Fi is, however, two to one: the invariance group of these points is bigger, Z 2 x2 C G. 

• The hyperbola 

4-t^ = 1 (44) 
ai ai - ai 

in the g 2 = plane passing through F 2 and A and having foci F 3 (the edge A 2 = o\ = A3 
in /MO)). 

The above points and curves play a special role in the definition of the elliptic coordinates 
and are also "critical loci" of the dynamics. 

3.2 Generic Kinks 

The generic kinks of the N = 3 system are the trajectories given by the solutions of (j35|)- (|36|) 
for non-zero finite values of C 2 and C 3 . The solutions of the implicit equations ([55|)-(|3(J) 
cannot be graphically represented by means of the built-in functions of Mathematica. We 
use a numerical algorithm implemented in Mathematica to obtain the graphic portrait of 
the trajectories. The algorithm allows us to calculate an arbitrary number of points on the 
orbit. These points joined by straight segments provide a visualization of the trajectory. 
There is a special step and an iteration of routine steps in the procedure, which is based on 
the Newton-Raphson method. 

First step. Identification of two points on the trajectory. 

For given values of C 2 , C3, cr 2 , cr 3 and a choice of signs, we set the first variable to the 
"point" Ai = Ai. (H)-H) becomes a system of two equations in two unknowns that can be 



19 



solved by the Newton- Raphson method with starting values (A 2 , A3). The outcome is a point 
Pi = (Ai, A2, A3) on the trajectory. We repeat this operation starting from Ai = A~i + e = \[ 
to find a second point P 2 = (A' 1; A' 2 , A 3 ) on the orbit. 

Ai, A' 1; \9, and A3 are chosen at random; good convergence is attained if these points 
belong to the middle zones of the variation ranges of Ai, A 2 and A3 or, at least, they are far 
away from the singularities on the faces of P^O) 
Successive steps. 

Pi and P2 provide an approximation of the curve by the secant line joining them. For 
some small 5 G R + , we choose P® = Pi + <5(P 2 — -Pi) as the starting value of the Newton- 
Raphson procedure applied to the solution of equations (|35|) and (|36|); we thus obtain the 
point P3 on the curve. P 2 and P 3 lead to guesstimate by the same token another value P° 
that produce the next point P4 on the orbit and now the iteration is obvious. Replacing S 
by —5 we travel along the opposite sense on the trajectory. The algorithm stops when one of 
the three variables Ai, A 2 , A3 reaches its extreme value; it is applied independently on each 
stage, determined by the signs of 7i a and the global trajectory is obtained by the demand for 
continuity. 

We now describe the portrait of these orbits. Having fixed 7 2 and 73, the corresponding 
kink trajectory is a non-plane curve in the interior of P^O) that starts from the vacuum point 
D, reaches the top face BCFiO and hits the edge AF 2 . It then goes to the edge F1F3, back 
again to the top face, hits the edge AF 2 a second time, the top face a third time and ends 
at D: see Figure 2 and Figure 3. Varying 7 2 and 73 in the range of finite real numbers, other 
similar trajectories are obtained that hit the edges AF 2 and F1F3 at different points. Given a 
sense of time there therefore exists a two-parameter family of kink trajectories in one-to-one 
correspondence with the points in the interior of AF 2 and F1F3. It should be mentioned 
that a whole congruence of trajectories parametrized by the interior of F1F3 converges at 
one single point in the interior of AF 2 and viceversa. 




Figure 2: A generic kink drawn both in D ?J and in ^3(0). Observe in D 3 that the generic kink is 
a heteroclinic trajectory 

The translation of a generic kink trajectory to Cartesian coordinates is a delicate matter; 
due to the non-uniqueness of the mapping implied by the change of coordinates special care 
is necessary in the analysis of the trajectory near the special conies P5[)-(PD where several 
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Figure 3: Several generic kinks, all of them intersecting once with the edge F1F3 and twice with 
the edge AF 2 of P 3 (0). 



options are a priori possible. Since it requires continuity and derivability to the trajectories 
in the interior of the ellipsoid ([40D, the interior of D 3 , the behaviour of the curve is mostly 
fixed. Two choices, the specification of D= v~ as the starting point and the location of the 
intersection of the kink trajectory with the q\ = plane in the quadrant characterized by 
Q2 > 0, 93 > 0, completely fix the itinerary. 

There is a crossroad where the "particle" touches the q± > branch of the hyperbola fl4*4]) 
and turns back towards the ellipse ([43]). There, the movement enters the q% < half-space 
and the kink trajectory reaches the other branch, q± < 0, of the critical hyperbola after 
crossing the g 2 = plane. At this stage, the particle makes its way for a third crossing of 
the q2 = plane and, finally, the journey ends at D= v + . This kind of kink trajectory is 
therefore heteroclinic: it starts and ends at different unstable points, so that: 

Ql = \ II ^ = ±1 ' <# = <# = (45) 
We call these topological kinks TK3 because they have three non-null components: 
qi {r) = fcix) ± 0, q 2 {r) = 2 (x) ^ 0, q 3 (r) = fa(x) ± 

It should be noted that a unique, apparently non-derivable, kink trajectory in elliptic 
coordinates corresponds to eight derivable trajectories in Cartesian coordinates: the choices 
of v~ or v + and q^ < or g 3 > 0, g 2 < or q 2 > as the starting point and initial quadrant 
give the eight possibilities. 

The energy of a three-component topological kink is the action of the trajectory times 

3 1— 

■^rm and hence computable from formula ( p9|) for the = 3 case: 



A 2 x/2 



E 



TK3 



\1d\1 



y/1 
4 2 

3 + 3 



Ai 



+ 



2A 2 rfA 2 



+ 



1 3A,rfA, 



0-3(3 - of) + o 2 (3 - o. 
It is independent of 7 2 and 73 and hence the same for every kink in the TK3 family. 



(46) 
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3.3 Enveloping Kinks 



There is another family of N = 3 kinks living on the surface M 3 = {(Ai, A 2 , A 3 ) / X± — 0}, the 
unique face of ^(O) where the elliptic coordinates are not singular. In M 3 , the Hamiltonian 
becomes: 



H 



E 



2A(X a ) 2 

-7T„ 



l Ag(A Q -a 2 2 )(A Q -a g / 
2 



(47) 



a=2K A'(A a ) a 2 A'(Aa) J 

and therefore there is a two-dimensional system hidden inside the A = 3 model which is 
Hamilton- Jacobi separable. The orbit equations, 



C 



VT^A, 



<7 2 



VI - A 2 + 



cr 2 



Vl — A3 - 0"2 



0"2 



o-3sign(7r 2 ) 



o- 3 sign(7r 3 ) 



VI - A 2 + (T 3 



\/l - A 2 - cr 3 



Vl - A 3 + 0-3 



Vi - A 3 - 0-3 



CT 2 sign(7r2) 



o- 2 sign(7r3) 



(46 



are parametrized by only one real constant 7 2 (C = e " 20 " 3 ^ ^D'w). 



Figure 4: The NTK3 family 

The Mathematica plot of these solutions is shown in Figure 4. Having fixed 72, the 
corresponding kink trajectory is a plane curve in M 3 that starts from the vacuum point 
D, reaches the top edge BC, goes to the umbilicus A and then back to the edge BC, to 
end finally in the vacuum point D. The value of 7 2 determines the points in BC where the 
trajectory bounces back and thus the one-parameter family of this kind of kink trajectories 
is in one-to-one correspondence with the points in the interior of BC. 

In Cartesian coordinates the enveloping kinks are trajectories that unfold on the ellipsoid 

^2 ^3 

The starting point is either D= v + or D= v~ and the trajectories also end in either 
D= v + or D= v~ . Associated with "homoclinic" trajectories, the corresponding kinks are 
"non-topological" : Q\ = Q\ = Q% = 0. The three Cartesian components q a differ from 
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zero and the appropriate name for this kind of solitary wave is a non-topological kink of 
three components, NTK3 for short. Every NTK3 trajectory on its way from D= iP 1 to 
D= -iP 1 crosses the umbilicus point of the ellipsoid. Note that, again, eight trajectories in 
the Cartesian space R 3 correspond to one trajectory in ^(O): the particle has the freedom 
to choose the points v + or v~ as base points of the curve. Having fixed one of them, the 
trajectory may develop in the half-ellipsoids determined in ( P5| ) by q 3 < or q 3 > and, 
finally, there are two travelling senses in each orbit. 

Also, the energy of a three-component non-topological kink is essentially the action of 
the NTK3 trajectory: 

AV2 fl \ 2 d\ 2 ri 2\ 3 d\ 3 J, n + uq] 

according to the Hamilton- Jacobi theory. 



3.4 Embedded Kinks 

Three-component topological and non-topological kinks arise as genuine solitary waves in 
the N = 3 model. Restriction to the q 3 = and/or q 2 = planes shows that the N = 2 
system is included twice, once in each plane, in the N = 3 model. Therefore, all the solitary 
waves of the N = 2 model are embedded twice as kinks of the larger N = 3 system. The 
embedded kinks live on the q 2 = and q 3 = planes, i.e. the faces of ^(O) where the elliptic 
coordinate system is singular. 



I. Embedded kinks in the q 2 = plane 

Both A 2 = o\ and A 3 = of give q 3 = 0, see (p8|), and hence this coordinate plane in R 3 is 
the union of the two faces, A2 = o\ and A3 = erf, of ^3(0). Therefore, in 

M 2a3 = {(A 1; A 2 , A 3 )/A 3 = al) U {(Ax, A 2 , A 3 )/A 2 = a 2 2 } = U 

we expect to find all the kinks of the N = 2 case. 

• In M 2 ^, A3 = of, we are in the face of -P^O) such that < Ai < of < A 2 < of, and 

q\ = -4(1 - A0(1 - A 2 ) , q 2 3 = 4fe 2 - Ai)^ 2 - A 2 ) 



The Hamiltonian also reduces to the N = 2 Hamiltonian 



2A(Xi) 2 1 / 2 2 . 2 

-n = , 7- . 7T-, , ; - ; A, Ai — cr 9 Ai — (Jo 

A'(Ai) 1 2A'(Ai) V 1V 2A 3 



2A(A 2 ) 2 
A'(A 2 ) 2 2A 



-^(^(^-^(A.-cr 2 )) 
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and the Hamilton- Jacobi method prescribes the equation 





Ai 


- 0"3 




Vi 


-Ax 


+ 1 


Vi- 


Ai 


+ 0"3 




Vi 


-Ai 


- 1 


Vi- 


A 2 


- 0"3 




Vi 


-A 2 


+ 1 


Vi- 


A 2 


+ ^3 




Vi 


-A 2 


- 1 



0-3 \ sign(vri) 



0-3 \ sign(7r 2 ) 



(50) 



as ruling the portion of the trajectories at this face, bounded by the edges AD, AF 2 , 
F 2 F 3 and F 3 D (see Figure 5). 





Figure 5: Several NTK2 a:j kink trajectories of the N = 2 system embedded in the g 2 = plane 
• In Mf , A 2 = o\ and the face in the boundary of £3(0) is < Ai < of < a\ < A3 < 1. 



The trajectory equations at this face are: 

Vl - Ai - 0-3 



32 (7 3 ct|72 



VI - Ai + a 3 



VI - A3 - 0-3 




0-3 \ sign(TTi) 



0-3 \ sign(7r 3 ) 



VI - A3 + CT 3 



and the boundary is formed by the edges AF 2 , F 2 0, OB and BA. 



(51) 



For finite values of 72, some kink trajectories given by ((50|)-([5l|) are depicted in Figure 
5. It may be observed that the trajectory starts at D and then runs through the face 
until the edge AF 2 . From this point, the particle enters the face (here the path is not 
derivable), reaches the BO edge and comes back to the AF 2 edge. This is the second point of 
non-differentiability re-entering the trajectory the M\ a face. All the trajectories then meet 
at the vertex F 3 and come back in a symmetric way to end in the D point. In Cartesian 
coordinates, these kink trajectories start and end in either D=v + or D= v~ , do not leave 
the g 2 = plane, and cross either the focus (qi = — cr 3 , q 3 = 0) or (qi = <t 3 , g 3 = 0). We 
therefore call them NTK2 CT3 because they are two-component non-topological kinks, merely 
the family of NTK2 kinks of the N = 2 model, embedded this way within the manifold of 
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kinks of the N = 3 system. There are four trajectories of this kind inside the ellipsoid fl4"5] ) 
in R 3 per trajectory in the boundary of -Ps(O): there is freedom to choose v + or v~ and 
the sense of travel in each orbit. The NTK^g kinks are fixed points of the Z2 sub-group of 
G = 3 generated by q 2 — > — c/2, that, however, does not leave invariant the TK3 and the 
NTK3 trajectories . The energy of these solutions is: 



\ 2 V2 



E 



1 2A,rfA, 



77T 



NTK2 




(52) 



There is a limiting case to this family of kinks: a trajectory along the DA and AB edges 
and back to D through the same way. Elliptic coordinates are even more singular on the 
edges, but the dynamical system reduces to a one-dimensional Hamiltonian system which 
can be integrated analytically. We have a two-step trajectory: 

At the DA edge, Ai = and A3 = of, the canonical equations (after use of the first 
integral) reduce to: 

d\ 2 



with the solution 

. TK2 CT3 

Ai 3 1 



T 



dr 



aJ K2ct3 1 



±2(A 2 -^)Vl-A 



T 



1 — erf tanh 2 (<73T) 



an 



(53) 



(54) 



for r G (— 00, arctanh ^] U arctanh 00). The second step occurs on the AB edge, 



<?3 



&3 



where, again, the canonical equations reduce to a single differential equation: if Ai = and 

A, 



^2 



- °"2> 

has the solution 

TK2 



dr 



aJ K2 -, 



±2(A 3 -^)Vl-A 



T 



a. 



2' 



aI K2ct3 1 



r) = 1 — af tanh 2 (cr 3 r) 



(55) 



(56) 



for t G 
TK2, 



Ai 3 (^)=0, _ 

— arctanh — , — arctanh — . The corresponding kinks in Cartesian coordinates are 

_ C3 ff3' (73 0-3 J ^ 

and TK2* the four two-component topological kinks of the N 



2 model: 



/ TK2 CT , . 
/ ?1 ( r ) 












TK2 CT3 , , 









tanh((T3T) 


±cr 3 sech(cr3r) 



V 



,TK2 
,TK2 

h 

,TK2 
^3 



X,t) \ 
£,t) / 



tanh(<T3a;) 
±1 ' 

±a 3 sech(cr 3 a;) 



Thus, the enveloping kinks of the N = 2 model are also embedded in the N 
The energy for these solutions and their anti-kinks is: 



\ 2 V2 



E 



rrr 



TK2 D 



\2d\2 



+ 



1 A 3 g?A 3 
1 x/T^Aa 



2(7, 1 



On 



(57) 
3 system. 



(5f 



In the c/2 = plane there is still one trajectory that is even more singular: it is a three 
step path running on the edges DF 3 , F 3 F 2 , F 2 and back to D through the same way. The 
canonical equations and its solutions in the three steps are: 
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1. A2 = o\ and A3 = of. 
d\ x 



dr 



= ±2Aiv/l — Ai, r G (— 00, — arctanha 3 ] U [arctanh <r 3 , 00) 



A™ 1 (r) = l-tanh 2 r ! Aj K1 (r) = <x 3 2 , Aj K1 (r) = a 2 
2 „„ j \ . _ ^2 



2. Ai = <r 3 and A3 = o" 2 
c?A 2 



±2A 2 y / l — A 2 , r G [— arctanho" 3 , arctanho" 2 ] U [arctanhcr 2 , arctanh<7 3 ] 



A* K1 (r) = al A 2 iK1 (r) = 1 - tanh 2 r, A 3 iK1 (r) = a, 



3. Ai = o\ and A 2 = a\. 

d\ 3 



dr 



= ±2A 3 y / l — A3, r G [arctanh<7 2 , arctanh(X 2 ] 



A^W^s 2 . A 2 rK1 (r) = a 2 2 , Aj K1 (r) = l-tanh 2 r 
Only one Cartesian component is different from zero: 





1 ±tanhr \ 


( <t>J K \x,t) \ 




1 ± tanh x 


\ 







$ K1 (x,t) 


- 







\ vJ k1 (t) J 


v ; 


Wl K1 0M) J 




V 


/ 



(59) 



and hence the one-component topological kink of the N — 1 model is embedded first in the 
manifold of kinks of the N = 2 model, and then in the iV = 3 system. There are two kinks of 
this kind in Cartesian coordinates which are mapped in a unique trajectory in the boundary 
of -P 3 (0). The TK1 trajectories are fixed points of the Z x2 sub-group of G generated by 
<?2 — > -<?2 and q 3 — > -q 3 . The energy is 



\ 2 V2 



TKl 



\\d\i 



Embedded Kinks in the q 3 = plane 



f I \ 2 d\ 2 



+ 



§ VI- A 2 ia| VI - A 



1 X 3 d\ 3 _ 4 
* = 3 



(60) 



The DF 3 , F 3 F 2 and F 2 edges form the intersection of the g 2 = and q 3 = planes. 
Therefore, the TKl kinks also live in the q 3 = plane. There is another maximally singular 
trajectory living on the "edge" in the q 3 = plane: 

At the DC edge, A 2 = a\ and Ai = 0, the canonical equations for the finite action 
trajectories are: 

L2(A 3 - d*)Jl^ (61) 



dr 



The path 



A! K ^(r)=0, \ 2 LK ^(r)=al (r) = 1 - aj tanh 2 (a 2 r) 



(62) 
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solves ( |6"I| ) and runs when r goes from — oo to +00 from D to D passing through the vertex 
C at t = 0. In Cartesian coordinates we recover the four two-component topological kinks 
of the N = 2 model, now embedded in the q 3 = plane: 



/ TK2 CT „ 
/ Ql 


(t) 






TK2 CT2 


(r) 






TK2 CT2 








V 93 


(r) 





tanh(cr2r) 
±af sech(cr2r) 




V 



(x,t) \ 
(x,t) 



,TK2 D 
>1 

,TK2 
>2 



± 



tanh(cT2a;) 
±d"2 sech(cr2x) 




(63) 

These are heteroclinic trajectories that produce the TK2 CT2 and TK2* 2 topological kinks. 
The energy is: 



E 



TK2„ 



1 AsdAs 



2cro 1 1 — -zr 



(64) 




Figure 6: NTK2 CT2 kink trajectories of the N = 2 system embedded in the (73 = plane 

Of course,the full manifold of kinks of the N = 2 model is embedded in the q% = plane 
: the set of kinks of the N = 3 system is completed by the two-component non-topological 
kinks living in the q 3 = plane, see Figure 6. The q 3 = plane is mapped to the union of 
two faces in the boundary of £3(0): 

M 2a2 = {(A 1; A 2 , A 3 )/A! = a 2 } U {(A l5 A 2 , A 3 )/A 2 = a 3 2 } = U 

• 1. In Ai = a 2 implies q 3 = 0. Therefore: 

a 2 

12 = —(^2-^2)(^2- A 3) 
0"2 

is a well defined change of coordinates in the range a\ < A2 < o\ < A3 < 1. In this 
region, the interior of the ellipse (f43|), the trajectories providing kinks are given by the 



27 



equations: 




(X2\ sign(7T2) 



a 2 \ sign(7r 3 ) 



(65) 



2. In , A 2 = 0-3 also implies g 3 = 0. In the range < Ai < a 3 < o\ < A 3 < 1 the 
change of coordinates is defined as 



9? 



9? 



-a 2 2 



:i-Ai)(l-A 3 ) 
^2 ~ Ai)(af - A 2 



The kink trajectories satisfy the equations: 



D 2cr 2 o-|72 



\/l - Ai - CT 2 



VI - Ai + cr 2 

VI - A 3 - 0"2 



VI - A 3 + cr 2 




0"2 \ signal) 



a 2 \ sign(7r 3 ) 



(66) 



The features of this kind of kinks are identical to the characteristics of the two-component 
non-topological kinks that exist in the = plane. The only difference is that they have 



support in the faces M\ 



and M? 



instead of Ml 

CT 2 *<T3 



and M? and we therefore call them 



NTK20-2. They meet at the vertex F 2 , and therefore at the foci {q\ = ±cr 2 ,0, 0) in R 3 ; see 
Figure 6. The energy is: 



A 2 V2 



E 



rrr 



NTK2 a 




(67) 



In sum: the manifold of kinks of the N = 2 model is embedded twice in the N = 3 
system, once in the q2 = plane and other in the g 3 = plane. They are sewn togehter by 
the common TK1, embedded from the N — 1 model . The embedded kinks fill the gaps left 
by the TK3 families of kinks in the interior of the ellipsoid fl45|) and also develop through 
the curves left by the NTK3 families on the boundary of D 3 . D 3 is thus a "totally" geodesic 
manifold with respect to the separatrices between bounded and unbounded motion in the 
N = 3 dynamical system. The NTK3 family form the envelop of the separatrices and the 
NTK3 kinks are themselves separatrices. 
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Figure 7: Plot of the "singular" topological kinks TK2 CT2 (solid line), TK2 CT3 (broken line) and 
TK1 (dash-dotted line) in both Cartesian and Elliptic coordinates. 

4 Further Comments 

We now infer the general structure of the kink manifold of the linear 0(iV)-sigma model 
from the pattern shown by the 0(2)- and 0(3)-sigma models. We can safely state that all 
the kink trajectories live in the sub-manifold D N C H N determined by the inequality: 

g? + 4 + - + 4< 1 

There are three categories: 
1. Generic Kinks 

A. There exists a family of generic kinks parametrized by iV — 1 real constants that 
live in the interior of D N . The intersection loci of the generic kinks are the singular 
quadrics: 

# + + =1, 9n = 0; Xi = \2 = a% 

U N u N "2 U N u N-1 

+ ••• -iA- =1, Qn-i = 0; A 2 = A 3 = a%_, 



'N-1 U N-1 u 2 "N-1 "N-2 



= 1 q 2 = 0; Aat_i = \ N = <j% 

"2 "3 "2 "N 2 

B. Generic kinks are non-topological- hence NTKN- if N even, and topological- hence 
TKN- if N is odd. 

2. Enveloping kinks. 

A. The restriction of the dynamical system to the boundary dD N of D N , the hyper- 
ellipsoid: 

g 1 2 + | + ... + 4 = l 

ai aft 
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provides a family of enveloping kinks parametrized by N — 2 real constants. Re- 
calling that in elliptic coordinates dD N is characterized by the equation Ai = 0, 
the intersection loci of this congruence are the umbilical sub-manifolds: 

Ai = , \ a = a 2 N ^ a+l = \ a+ i , Va = 2,3, ...,N - 1 

of dimension iV — 3 of the hyper-ellipsoid dD N . 

B. Enveloping kinks are topological, hence TKN, if N is even, and non-topological, 
hence NTKN, if N is odd. 

3. Embedded Kinks. 

On the N — l R^" 1 sub-manifolds determined by the conditions q a = 0, if a = 2 or 3 
or ... or N, the dynamical system reduces to the mechanical system that arises in the 
linear 0(N— l)-sigma model. Thus, the kink manifold of the N—l case is included N—l 
times in the 0(iV)-model, filling the holes left in the interior of D N by the generic kinks, 
and also covering in dD N the sub-spaces which are not covereded by the enveloping 
kinks. Each N—l kink sub-manifold is not, however, included (N — 1) x (N — 2) 
times in the 0(iV)-model because the R^" 2 sub-spaces are intersections of the N — l 
R^ 1 , defined above. The N — l kink manifolds are not separated but sewn togehter 
through the N — 2 kink sub-manifolds. This is a iterative process in such a way that 

the kink manifold of the 0(N — r)-sigma model is included ^ ^ ^ ^ ^ times in the 

kink manifold of the 0(N) system. 

One can ask what happens if a continuous sub-group 0(r) of O(N) survives as symme- 
try group of the system. This happens if the deformation is chosen in such a way that 
< a\ — of = • • • = a 2 < a 2 +1 < . . . < a 2 N < 1. In this case we obtain a sub- manifold 
of kinks from 0(r) rotations around the q\ axis of the kink manifold of the N = 2 
system that lives in the q\ : q<i plane. The remaining kinks correspond to the solitary 
waves of the N = r — 1 system defined in the orthogonal R Ar - r + 1 sub-space. Also 
the deformations where 1 < of +1 < ■ ■ ■ < <j 2 n are easy to understand. Finite action 
trajectories spread out in the domain in H N bounded by the hyper-hyperboloid: 

2 , , , Qr ll+l _ _ In _ 

l/l-r_2~r'''~r_9 1-2 i ' ' ' |-2I — 
Oi 2 r Wr+l\ Kl 

The kink manifold of this system is the kink manifold of the N = r model defined in 
the sub-space R r C such that q r+ \ — ■ ■ ■ — qN — 0. 

Finally we consider a mild deformation of our model by introducing asymmetries in 
the non-harmonic terms of the potential energy and also adapting the quadratic terms 
in a suitable manner: 

V = i(0? + (l+£ 2 )02 + ... + (l + £ JV )0^-l) 2 + 



+ y0 2 + - + ~y<Pn + J<P2 + - + -j<Pn 
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The new non-dimensional constants e a and 5 a are defined in terms of the old o~ a 's 
through: 

1+e „ = _>|_J>, S a = -^\ +a '\ . = 2,3,..., AT 

Among the kinks of the deformed linear 0(iV)-sigma model only the following survive 
as solitary wave solutions of this perturbed system: 

(a) The TK1. 

</>i = tanhx , (f>2 = ... = 4>n = 



(b) All the TK2 kinks. On the ellipse, 



2 _|_ 1 ~l~ £ a 2 



the TK2a a and TK2*cr a configurations 



1 - CT 2 

0i = tanha a a; , a = ±W 21 secha a x , 0^ = 0, V b ^ a,b ^ 1 

V 1 +£a 

are solutions of the field equations. The amazing fact is that in this deformation 
of the 0(iV)-linear sigma model, discussed by Bazeia et al. if iV = 2 [T2j, the 
energy of all these kinks is the same: 

Am 3 

E TK i = E TK2 a 2 — ■■■ — E T K2a N = ^/^^2~ 

On one hand, we have a deformation of the linear 0(iV)-sigma model that exhibits 
a complex variety of kinks; on the other hand, another deformation of the O(N)- 
model rejects almost every kink but the simplest ones as solutions, and all of the 
surviving kinks are degenerated in energy. 



5 Outlook 

The developments disclosed in this paper suggest a general strategy in the search for kinks 
in two space-time dimensional field theories. When the fields have N components assembled 
in a vector representation of the O(N) group, we focus on systems with symmetry breaking 
to a discrete sub-group of O(N) which has more than one element. If the dynamical system 
that determines the localized static solutions is completely integrable, all the solitary waves 
can be found, at least in principle. Particularly interesting is the situation where the N — 1 
invariants in involution with the mechanical energy act non-trivially on the manifold of 
localized solutions and the orbit is a continuous space. One can then perturb such a system, 
loosing in the perturbation many of the solitary wave solutions: only few of the localized 
static solutions survive as kinks of the perturbed (more realistic) model. 

We finally list several interesting questions that will be postponed for future research: 

• study of the structure of the kink manifold of the deformed linear 0(iV)-sigHia model 
as a moduli space seems to be worthwhile. 
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• A detailed analysis of the sum rules between the energies of the different kinds of kinks 
is necessary to fix the structure mentioned above. 

• A treatment a la Bogomolny is also possible. This allows for a supersymmetric exten- 
sion of the model in such a way that the kinks become BPS states. 

• The difficult problem remains of determining the stability of the different kinds of kinks. 
Application of the Morse index theorem helps in finding the stability properties, which 
in turn provide information about the quantization of these topological defects. 
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Appendix: Elliptic coordinates 

Given any set of N real positive numbers such that < r\ < r 2 < . . . < rjv, let us consider 
the equation: 

N 2 

= ! (68) 

a=l ' a A 
N g 2 

The left-hand member Q\(q) = \^ — of this equation can be understood either as a 

function of R , for fixed A G C, or as a function of the complex variable A, for fixed q G R w . 
From fl68|) one immediately deduces: 

N , fi(A-Aa) 

1 " Qx(q) = 1 + £ = ^ = ( 69 ) 

0=1 a II(A-r.) 

a=l 

Therefore, the N roots A a of the polynomial in the numerator of 1 — Q\(q), a rational 
function of A, are the roots of equation (|68|) . The roots A a are also real numbers and 1 — Q\{q) 
is a rational function such that Ai < r\ < A 2 < . . . < r^-i < Ajv < r N , see Figure 8. To 
prove this point one needs to study 1 — Q\(q) along the A-real axis, near the poles A = r a , 
using Bolzano's theorem. 

Definition. The elliptic coordinates of the point q = • ■ ■ , <7jv) £ R w are the roots 
\ E = (A 1; ...,X N )e P N (oo) of Q\(q) = 1. 

P N (oo) C H N is the open sub-space of given by: — oo < Ai < r\, r\ < \ 2 < r 2 , . . ., 
r^v-i < Aat < r N . The solution of ( |5ED for A = Ai =constant G (— oo,ri) is, geometrically, a 
quadric surface, a hyperellipsoid, when q varies in R^. The family of quadrics obtained by 
taking A = A a = constant G (r a _ 1; r a ), a > 2, corresponds to a family of hyper-hyperboloids 
of every possible signature in R^, if Q\(q) is considered as a function of q. 
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3 „2 

/•„ - A 



Figure 8: Plot of the function y = ^ 2—^- and the function y = 1, see (|68[), for fixed values of 

<? a and r a . 



a=l 



It is convenient to denote the products in (^) as: 

N N 

A(A) = n(A-A a ), A(A) = n(A-r.) 



so that 



iV 2 
1 + E '" 



N 
a=l 



^ A - r ^ 

11 (A — r c 

a=l 



A(A) 
A(A) 



(70) 



An explicit formula for defining q\ as a function of the A's , Va, is obtained by applying 
the residue theorem to both members of the equation (ffOp: 



Res(l-Q A (q))| 



A=r a 



(71) 



A(ra) -, A'(r a ) = dA(A) 



A=r a 



A'fr„ 



dA 



A=r„ 



Therefore: 



92 



A(r c 



v 

II( r a-A6) 

b=l 

A'(r a ) ~ * 

11 ( r a-n) 

b=l,b^a 



(72) 



and we see that the transformation (gx, . . . , q^) — ► (Ai, . . . , Ajv) is 2 W to 1. 

Inverting ( [72] ) to express A a as a function of the g's, Va, requires one to solve an algebraic 
equation in A a with powers up to A^. This is easy for N = 2, possible, but very difficult 
for N = 3,4 using Cardano's formulas, and impossible if N > 5. For this reason, another 
derivation of (^) is useful, which in passing allows one to show identities between Cartesian 
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and elliptic coordinates that make practical computations possible. To do this, notice that 
(|70| ) implies: 

N N N 

A(A) = n(A-r a ) + E^ n ( A "^) ( 73 ) 

a=l a=l b=l,b=£a 

Setting A = r c in (73) one immediately derives (|72|) . More important, expanding the two 



members of (|73|) in a power series in A , we obtain 

TV N \ I N N N 



/ N N \ / N N N \ 

v v + a^ -E^ + Eq+^ 2 EE^-E^ E n )+... 

\ a=l a=l / \ b<a a=l b=l,b^a J 

/ N \ ( N \ N 

= * N - E M ^ + E E a^ 2 + . . . + (-if n a. 

\a=l / \ b<a J a=l 

Equalizing the coefficients of the terms with the same power of A in the last equation we 
have iV non trivial identities. We shall use the equalities between the coefficients of A^ -1 
and A^" 2 : 

N N N 

= (74) 

a=l a=l a=l 

N N N N N 

E E X * X b = E r «<?a - E <?a E r b + E E r b r a (75) 
b<a a=l a=l 6=1 b<a 

Another important tool using elliptic coordinates is the Jacobi lemma: 
Lemma. The expression 



N a s 

1 'a 



E 



(a a - ai)(a a - a 2 ) ■ ■ ■ ...(a a - a N ) 

where (a%, . . . ,«tv) are A" real numbers such that «i < a 2 < ■ ■ ■ < (%n is equal to if 

s < N- 2 and 1 for s = N- 1. 

Proof: 

Consider the function 

fs(z) = ^ 

T[(z-a a ) 

a=l 

which has A" poles in the complex plane at z = a a , Va, and another pole at infinity. If 7 is 
a closed curve which is the boundary of a region D of the complex plane containing all the 
finite poles of f s (z), the residue theorem tells us that: 

1 / Jt If 

— & f s {z)dz = £ Res(/ S )(a a ) = -Res(/ s )(oo) = - — f f s (z)dz 
2m J-y ~y 2m J-j 



Also, 



N N 



ERes(/ s )(a a ) = E 
0=1 

Res(/,)(oo) 



t^i " «^=i (ot a - ai)(a a - a 2 ) . . . . . . (a a - a N ) 

0, Vs < A^ — 1 
-1, s = N-l 
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and the lemma is proved. 

Applying this result to the choice a a — A , Vo = 1, . . . , N, we obtain new identities 

N \s 
N \N-1 

Sa^T' (77) 

because A'(A a ) = (A — Ai)(A — A 2 ) . . . SV. . . . (A a — Ajy). Alternatively if we take a a = r a , 
the lemma implies that: 

N s N N-l 

y-rr L T = 0,Vs< N-l; Y-2— = l (78) 
An important identity obtained from the lemma is: 

N 2 

E fr Y%Z X) =°^b,c (79) 
The normal vectors to the family of quadrics (|68|) 

N 2 

QM = Y. 1 



o=l r a - A 



at the point q = (g 1; . . . , q N ), are 

<?i 9jv 



n(A) = (ni(A), . . . ,njv(A)) 
Observe that (|79|) implies 



r\ — A r^v — A 



v 

^n a (Afc)n a (A c ) =0, Va,6 

0=1 

Therefore, all the quadrics are orthogonal with each other and the elliptic coordinates form 
an orthogonal system. The standard Euclidean metric in Cartesian coordinates can be 
expressed in elliptic coordinates in the form: 

N N N 

ds 2 = J2dql = J2J2 9ab(^E)d\ a d\ b 

a=l a=l 6=1 

Derivation of the two members of equation (f72|) leads to: 



2dq a t ^ N d\b 



q a b=1 r a - h 
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and, using the Jacobi Lemma, 

N rl\ 2 

4drf = rfE b 



6=1 T 'a ~ \ 



Finally, we have: 



N 

„ N II (Aa"A 6 ) 

1 A (A a ) 1 6^a,ft=l w , , 

^ = -i^pg = -4-^ > = 0, Va ^ 6 

Il(Aa-r b ) 

6=1 

The kinetic energy of a natural dynamical system in elliptic coordinates 

1 N i N 

z a=l Z a=l 

In terms of the canonical momentum n a = J2^, T reads: 

n ^ N i N a(\ \ 



a=l 
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